Ketamine can produce oscillatory dynamics by engaging mechanisms dependent on the kinetics of NMDA receptors.

Ketamine is an NMDA-receptor antagonist that produces sedation, analgesia and dissociation at low doses and profound unconsciousness with antinociception at high doses. At high and low doses, ketamine can generate gamma oscillations (>25 Hz) in the electroencephalogram (EEG). The gamma oscillations are interrupted by slow-delta oscillations (0.1-4 Hz) at high doses. Ketamine's primary molecular targets and its oscillatory dynamics have been characterized. However, how the actions of ketamine at the subcellular level give rise to the oscillatory dynamics observed at the network level remains unknown. By developing a biophysical model of cortical circuits, we demonstrate how NMDA-receptor antagonism by ketamine can produce the oscillatory dynamics observed in human EEG recordings and non-human primate local field potential recordings. We have discovered how impaired NMDA-receptor kinetics can cause disinhibition in neuronal circuits and how a disinhibited interaction between NMDA-receptor-mediated excitation and GABA-receptor-mediated inhibition can produce gamma oscillations at high and low doses, and slow-delta oscillations at high doses. Our work uncovers general mechanisms for generating oscillatory brain dynamics that differs from ones previously reported, and provides important insights into ketamine's mechanisms of action as an anesthetic and as a therapy for treatment-resistant depression.

Multisensory gamma stimulation promotes glymphatic clearance of amyloid.

The glymphatic movement of fluid through the brain removes metabolic waste1-4. Noninvasive 40 Hz stimulation promotes 40 Hz neural activity in multiple brain regions and attenuates pathology in mouse models of Alzheimer's disease5-8. Here we show that multisensory gamma stimulation promotes the influx of cerebrospinal fluid and the efflux of interstitial fluid in the cortex of the 5XFAD mouse model of Alzheimer's disease. Influx of cerebrospinal fluid was associated with increased aquaporin-4 polarization along astrocytic endfeet and dilated meningeal lymphatic vessels. Inhibiting glymphatic clearance abolished the removal of amyloid by multisensory 40 Hz stimulation. Using chemogenetic manipulation and a genetically encoded sensor for neuropeptide signalling, we found that vasoactive intestinal peptide interneurons facilitate glymphatic clearance by regulating arterial pulsatility. Our findings establish novel mechanisms that recruit the glymphatic system to remove brain amyloid.

Rapid thalamocortical network switching mediated by cortical synchronization underlies propofol-induced EEG signatures: a biophysical model.

Propofol-mediated unconsciousness elicits strong alpha/low-beta and slow oscillations in the electroencephalogram (EEG) of patients. As anesthetic dose increases, the EEG signal changes in ways that give clues to the level of unconsciousness; the network mechanisms of these changes are only partially understood. Here, we construct a biophysical thalamocortical network involving brain stem influences that reproduces transitions in dynamics seen in the EEG involving the evolution of the power and frequency of alpha/low-beta and slow rhythm, as well as their interactions. Our model suggests that propofol engages thalamic spindle and cortical sleep mechanisms to elicit persistent alpha/low-beta and slow rhythms, respectively. The thalamocortical network fluctuates between two mutually exclusive states on the timescale of seconds. One state is characterized by continuous alpha/low-beta-frequency spiking in thalamus (C-state), whereas in the other, thalamic alpha spiking is interrupted by periods of co-occurring thalamic and cortical silence (I-state). In the I-state, alpha colocalizes to the peak of the slow oscillation; in the C-state, there is a variable relationship between an alpha/beta rhythm and the slow oscillation. The C-state predominates near loss of consciousness; with increasing dose, the proportion of time spent in the I-state increases, recapitulating EEG phenomenology. Cortical synchrony drives the switch to the I-state by changing the nature of the thalamocortical feedback. Brain stem influence on the strength of thalamocortical feedback mediates the amount of cortical synchrony. Our model implicates loss of low-beta, cortical synchrony, and coordinated thalamocortical silent periods as contributing to the unconscious state.NEW & NOTEWORTHY GABAergic anesthetics induce alpha/low-beta and slow oscillations in the EEG, which interact in dose-dependent ways. We constructed a thalamocortical model to investigate how these interdependent oscillations change with propofol dose. We find two dynamic states of thalamocortical coordination, which change on the timescale of seconds and dose-dependently mirror known changes in EEG. Thalamocortical feedback determines the oscillatory coupling and power seen in each state, and this is primarily driven by cortical synchrony and brain stem neuromodulation.

Rapid synaptic and gamma rhythm signature of mouse critical period plasticity.

Early-life experience enduringly sculpts thalamocortical (TC) axons and sensory processing. Here, we identify the very first synaptic targets that initiate critical period plasticity, heralded by altered cortical oscillations. Monocular deprivation (MD) acutely induced a transient (<3 h) peak in EEG γ-power (~40 Hz) specifically within the visual cortex, but only when the critical period was open (juvenile mice or adults after dark-rearing, Lynx1-deletion, or diazepam-rescued GAD65-deficiency). Rapid TC input loss onto parvalbumin-expressing (PV) inhibitory interneurons (but not onto nearby pyramidal cells) was observed within hours of MD in a TC slice preserving the visual pathway - again once critical periods opened. Computational TC modeling of the emergent γ-rhythm in response to MD delineated a cortical interneuronal gamma (ING) rhythm in networks of PV-cells bearing gap junctions at the start of the critical period. The ING rhythm effectively dissociated thalamic input from cortical spiking, leading to rapid loss of previously strong TC-to-PV connections through standard spike-timing-dependent plasticity rules. As a consequence, previously silent TC-to-PV connections could strengthen on a slower timescale, capturing the gradually increasing γ-frequency and eventual fade-out over time. Thus, ING enables cortical dynamics to transition from being dominated by the strongest TC input to one that senses the statistics of population TC input after MD. Taken together, our findings reveal the initial synaptic events underlying critical period plasticity and suggest that the fleeting ING accompanying a brief sensory perturbation may serve as a robust readout of TC network state with which to probe developmental trajectories.

Deep brain stimulation in the subthalamic nucleus for Parkinson's disease can restore dynamics of striatal networks.

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is highly effective in alleviating movement disability in patients with Parkinson's disease (PD). However, its therapeutic mechanism of action is unknown. The healthy striatum exhibits rich dynamics resulting from an interaction of beta, gamma, and theta oscillations. These rhythms are essential to selection and execution of motor programs, and their loss or exaggeration due to dopamine (DA) depletion in PD is a major source of behavioral deficits. Restoring the natural rhythms may then be instrumental in the therapeutic action of DBS. We develop a biophysical networked model of a BG pathway to study how abnormal beta oscillations can emerge throughout the BG in PD and how DBS can restore normal beta, gamma, and theta striatal rhythms. Our model incorporates STN projections to the striatum, long known but understudied, found to preferentially target fast-spiking interneurons (FSI). We find that DBS in STN can normalize striatal medium spiny neuron activity by recruiting FSI dynamics and restoring the inhibitory potency of FSIs observed in normal conditions. We also find that DBS allows the reexpression of gamma and theta rhythms, thought to be dependent on high DA levels and thus lost in PD, through cortical noise control. Our study highlights that DBS effects can go beyond regularizing BG output dynamics to restoring normal internal BG dynamics and the ability to regulate them. It also suggests how gamma and theta oscillations can be leveraged to supplement DBS treatment and enhance its effectiveness.

A biophysical model of striatal microcircuits suggests gamma and beta oscillations interleaved at delta/theta frequencies mediate periodicity in motor control.

Striatal oscillatory activity is associated with movement, reward, and decision-making, and observed in several interacting frequency bands. Local field potential recordings in rodent striatum show dopamine- and reward-dependent transitions between two states: a "spontaneous" state involving ß (∼15-30 Hz) and low γ (∼40-60 Hz), and a state involving θ (∼4-8 Hz) and high γ (∼60-100 Hz) in response to dopaminergic agonism and reward. The mechanisms underlying these rhythmic dynamics, their interactions, and their functional consequences are not well understood. In this paper, we propose a biophysical model of striatal microcircuits that comprehensively describes the generation and interaction of these rhythms, as well as their modulation by dopamine. Building on previous modeling and experimental work suggesting that striatal projection neurons (SPNs) are capable of generating ß oscillations, we show that networks of striatal fast-spiking interneurons (FSIs) are capable of generating δ/θ (ie, 2 to 6 Hz) and γ rhythms. Under simulated low dopaminergic tone our model FSI network produces low γ band oscillations, while under high dopaminergic tone the FSI network produces high γ band activity nested within a δ/θ oscillation. SPN networks produce ß rhythms in both conditions, but under high dopaminergic tone, this ß oscillation is interrupted by δ/θ-periodic bursts of γ-frequency FSI inhibition. Thus, in the high dopamine state, packets of FSI γ and SPN ß alternate at a δ/θ timescale. In addition to a mechanistic explanation for previously observed rhythmic interactions and transitions, our model suggests a hypothesis as to how the relationship between dopamine and rhythmicity impacts motor function. We hypothesize that high dopamine-induced periodic FSI γ-rhythmic inhibition enables switching between ß-rhythmic SPN cell assemblies representing the currently active motor program, and thus that dopamine facilitates movement in part by allowing for rapid, periodic shifts in motor program execution.

Biased competition in the absence of input bias revealed through corticostriatal computation.

Classical accounts of biased competition require an input bias to resolve the competition between neuronal ensembles driving downstream processing. However, flexible and reliable selection of behaviorally relevant ensembles can occur with unbiased stimulation: striatal D1 and D2 spiny projection neurons (SPNs) receive balanced cortical input, yet their activity determines the choice between GO and NO-GO pathways in the basal ganglia. We here present a corticostriatal model identifying three mechanisms that rely on physiological asymmetries to effect rate- and time-coded biased competition in the presence of balanced inputs. First, tonic input strength determines which one of the two SPN phenotypes exhibits a higher mean firing rate. Second, low-strength oscillatory inputs induce higher firing rate in D2 SPNs but higher coherence between D1 SPNs. Third, high-strength inputs oscillating at distinct frequencies can preferentially activate D1 or D2 SPN populations. Of these mechanisms, only the latter accommodates observed rhythmic activity supporting rule-based decision making in prefrontal cortex.

M-Current Expands the Range of Gamma Frequency Inputs to Which a Neuronal Target Entrains.

Theta (4-8 Hz) and gamma (30-80 Hz) rhythms in the brain are commonly associated with memory and learning (Kahana in J Neurosci 26:1669-1672, 2006; Quilichini et al. in J Neurosci 30:11128-11142, 2010). The precision of co-firing between neurons and incoming inputs is critical in these cognitive functions. We consider an inhibitory neuron model with M-current under forcing from gamma pulses and a sinusoidal current of theta frequency. The M-current has a long time constant (∼90 ms) and it has been shown to generate resonance at theta frequencies (Hutcheon and Yarom in Trends Neurosci 23:216-222, 2000; Hu et al. in J Physiol 545:783-805, 2002). We have found that this slow M-current contributes to the precise co-firing between the network and fast gamma pulses in the presence of a slow sinusoidal forcing. The M-current expands the phase-locking frequency range of the network, counteracts the slow theta forcing, and admits bistability in some parameter range. The effects of the M-current balancing the theta forcing are reduced if the sinusoidal current is faster than the theta frequency band. We characterize the dynamical mechanisms underlying the role of the M-current in enabling a network to be entrained to gamma frequency inputs using averaging methods, geometric singular perturbation theory, and bifurcation analysis.

Striatal cholinergic receptor activation causes a rapid, selective and state-dependent rise in cortico-striatal ß activity.

Cortico-basal ganglia-thalamic (CBT) ß oscillations (15-30 Hz) are elevated in Parkinson's disease and correlated with movement disability. To date, no experimental paradigm outside of loss of dopamine has been able to specifically elevate ß oscillations in the CBT loop. Here, we show that activation of striatal cholinergic receptors selectively increased ß oscillations in mouse striatum and motor cortex. In individuals showing simultaneous ß increases in both striatum and M1, ß partial directed coherence (PDC) increased from striatum to M1 (but not in the reverse direction). In individuals that did not show simultaneous ß increases, ß PDC increased from M1 to striatum (but not in the reverse direction), and M1 was characterized by persistent ß-high frequency oscillation phase-amplitude coupling. Finally, the direction of ß PDC distinguished between ß sub-bands. This suggests that (1) striatal cholinergic tone exerts state-dependent and frequency-selective control over CBT ß power and coordination; (2) ongoing rhythmic dynamics can determine whether elevated ß oscillations are expressed in striatum and M1; and (3) altered striatal cholinergic tone differentially modulates distinct ß sub-bands.

Thalamocortical control of propofol phase-amplitude coupling.

The anesthetic propofol elicits many different spectral properties on the EEG, including alpha oscillations (8-12 Hz), Slow Wave Oscillations (SWO, 0.1-1.5 Hz), and dose-dependent phase-amplitude coupling (PAC) between alpha and SWO. Propofol is known to increase GABAA inhibition and decrease H-current strength, but how it generates these rhythms and their interactions is still unknown. To investigate both generation of the alpha rhythm and its PAC to SWO, we simulate a Hodgkin-Huxley network model of a hyperpolarized thalamus and corticothalamic inputs. We find, for the first time, that the model thalamic network is capable of independently generating the sustained alpha seen in propofol, which may then be relayed to cortex and expressed on the EEG. This dose-dependent sustained alpha critically relies on propofol GABAA potentiation to alter the intrinsic spindling mechanisms of the thalamus. Furthermore, the H-current conductance and background excitation of these thalamic cells must be within specific ranges to exhibit any intrinsic oscillations, including sustained alpha. We also find that, under corticothalamic SWO UP and DOWN states, thalamocortical output can exhibit maximum alpha power at either the peak or trough of this SWO; this implies the thalamus may be the source of propofol-induced PAC. Hyperpolarization level is the main determinant of whether the thalamus exhibits trough-max PAC, which is associated with lower propofol dose, or peak-max PAC, associated with higher dose. These findings suggest: the thalamus generates a novel rhythm under GABAA potentiation such as under propofol, its hyperpolarization may determine whether a patient experiences trough-max or peak-max PAC, and the thalamus is a critical component of propofol-induced cortical spectral phenomena. Changes to the thalamus may be a critical part of how propofol accomplishes its effects, including unconsciousness.

Striatal cholinergic interneurons generate beta and gamma oscillations in the corticostriatal circuit and produce motor deficits.

Cortico-basal ganglia-thalamic (CBT) neural circuits are critical modulators of cognitive and motor function. When compromised, these circuits contribute to neurological and psychiatric disorders, such as Parkinson's disease (PD). In PD, motor deficits correlate with the emergence of exaggerated beta frequency (15-30 Hz) oscillations throughout the CBT network. However, little is known about how specific cell types within individual CBT brain regions support the generation, propagation, and interaction of oscillatory dynamics throughout the CBT circuit or how specific oscillatory dynamics are related to motor function. Here, we investigated the role of striatal cholinergic interneurons (SChIs) in generating beta and gamma oscillations in cortical-striatal circuits and in influencing movement behavior. We found that selective stimulation of SChIs via optogenetics in normal mice robustly and reversibly amplified beta and gamma oscillations that are supported by distinct mechanisms within striatal-cortical circuits. Whereas beta oscillations are supported robustly in the striatum and all layers of primary motor cortex (M1) through a muscarinic-receptor mediated mechanism, gamma oscillations are largely restricted to the striatum and the deeper layers of M1. Finally, SChI activation led to parkinsonian-like motor deficits in otherwise normal mice. These results highlight the important role of striatal cholinergic interneurons in supporting oscillations in the CBT network that are closely related to movement and parkinsonian motor symptoms.

Therapeutic mechanisms of high-frequency stimulation in Parkinson's disease and neural restoration via loop-based reinforcement.

High-frequency deep brain stimulation (HFS) is clinically recognized to treat parkinsonian movement disorders, but its mechanisms remain elusive. Current hypotheses suggest that the therapeutic merit of HFS stems from increasing the regularity of the firing patterns in the basal ganglia (BG). Although this is consistent with experiments in humans and animal models of Parkinsonism, it is unclear how the pattern regularization would originate from HFS. To address this question, we built a computational model of the cortico-BG-thalamo-cortical loop in normal and parkinsonian conditions. We simulated the effects of subthalamic deep brain stimulation both proximally to the stimulation site and distally through orthodromic and antidromic mechanisms for several stimulation frequencies (20-180 Hz) and, correspondingly, we studied the evolution of the firing patterns in the loop. The model closely reproduced experimental evidence for each structure in the loop and showed that neither the proximal effects nor the distal effects individually account for the observed pattern changes, whereas the combined impact of these effects increases with the stimulation frequency and becomes significant for HFS. Perturbations evoked proximally and distally propagate along the loop, rendezvous in the striatum, and, for HFS, positively overlap (reinforcement), thus causing larger poststimulus activation and more regular patterns in striatum. Reinforcement is maximal for the clinically relevant 130-Hz stimulation and restores a more normal activity in the nuclei downstream. These results suggest that reinforcement may be pivotal to achieve pattern regularization and restore the neural activity in the nuclei downstream and may stem from frequency-selective resonant properties of the loop.

Neurosystems: brain rhythms and cognitive processing.

Neuronal rhythms are ubiquitous features of brain dynamics, and are highly correlated with cognitive processing. However, the relationship between the physiological mechanisms producing these rhythms and the functions associated with the rhythms remains mysterious. This article investigates the contributions of rhythms to basic cognitive computations (such as filtering signals by coherence and/or frequency) and to major cognitive functions (such as attention and multi-modal coordination). We offer support to the premise that the physiology underlying brain rhythms plays an essential role in how these rhythms facilitate some cognitive operations.

Dynamical changes in neurological diseases and anesthesia.

Dynamics of neuronal networks can be altered in at least two ways: by changes in connectivity, that is, the physical architecture of the network, or changes in the amplitudes and kinetics of the intrinsic and synaptic currents within and between the elements making up a network. We argue that the latter changes are often overlooked as sources of alterations in network behavior when there are also structural (connectivity) abnormalities present; indeed, they may even give rise to the structural changes observed in these states. Here we look at two clinically relevant states (Parkinson's disease and schizophrenia) and argue that non-structural changes are important in the development of abnormal dynamics within the networks known to be relevant to each disorder. We also discuss anesthesia, since it is entirely acute, thus illustrating the potent effects of changes in synaptic and intrinsic membrane currents in the absence of structural alteration. In each of these, we focus on the role of changes in GABAergic function within microcircuits, stressing literature within the last few years.

Potential network mechanisms mediating electroencephalographic beta rhythm changes during propofol-induced paradoxical excitation.

Propofol, like most general anesthetic drugs, can induce both behavioral and electroencephalographic (EEG) manifestations of excitation, rather than sedation, at low doses. Neuronal excitation is unexpected in the presence of this GABA(A)-potentiating drug. We construct a series of network models to understand this paradox. Individual neurons have ion channel conductances with Hodgkin-Huxley-type formulations. Propofol increases the maximal conductance and time constant of decay of the synaptic GABA(A) current. Networks range in size from 2 to 230 neurons. Population output is measured as a function of pyramidal cell activity, with the electroencephalogram approximated by the sum of population AMPA activity between pyramidal cells. These model networks suggest propofol-induced paradoxical excitation may result from a membrane level interaction between the GABA(A) current and an intrinsic membrane slow potassium current (M-current). This membrane level interaction has consequences at the level of multicellular networks enabling a switch from baseline interneuron synchrony to propofol-induced interneuron antisynchrony. Large network models reproduce the clinical EEG changes characteristic of propofol-induced paradoxical excitation. The EEG changes coincide with the emergence of antisynchronous interneuron clusters in the model networks. Our findings suggest interneuron antisynchrony as a potential network mechanism underlying the generation of propofol-induced paradoxical excitation. As correlates of behavioral phenomenology, these networks may refine our understanding of the specific behavioral states associated with general anesthesia.